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Resumen 


La sequía es el fenómeno natural que afecta a más personas en el 
mundo, hasta la última consecuencia, la muerte de plantas, animales y 
seres humanos. Para evaluar los posibles daños que causa, se 
desarrollaron índices de sequía que facilitan su análisis cuantitativo para 
realizar una planeación y manejo de los recursos hídricos. Sin embargo, 
actualmente se observan cambios en el sistema climático global que 
provocan condiciones de sequía de mayor intensidad, duración y 
frecuencia. Por lo que se realizó la detección de cambios con el método 
de Rodionov en los índices de sequía SPEI (Índice de Precipitación y 
Evapotranspiración Estandarizado) y SPI (Índice de Precipitación 


Estandarizada) a escalas de 6 y 12 meses en el estado de Durango, 
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México. Se analizaron las tendencias de cambios en la media (RSI) y 
varianza (RSSI) de los índices en 49 estaciones del índice SPI de 6 
meses y en 45 del índice SPEI de seis meses. Los valores RSI 
identificaron el comienzo de los periodos de sequía, la mayoría de los 
valores máximos de RSI se presentaron en julio de 1992. Las tendencias 
de la frecuencia de valores RSI, disminuyó durante el periodo 2000- 
2004 y aumentó durante 2005-2008. La década de mayores cambios en 
la media fue 1990-1999. La mayor frecuencia de valores RSSI fue en el 
periodo 1985-1989, los cambios se mantuvieron constantes durante los 
últimos 20 años. La década de mayores cambios en la varianza fue de 
1980-1989. 


Palabras clave: Durango, sequía, SPI, SPEI, Rodionov, cambios de 


media, cambios de varianza. 


Abstract 


Drought is the natural phenomenon that affects the most people in the 
world, until the last consequence, the death of plants, animals, and 
humans. To evaluate the possible damages it causes, drought indices 
were developed that facilitate its quantitative analysis to carry out 
planning and management of water resources. However, there are 
currently changes in the global climate system, which cause drought 


conditions of greater intensity, duration, and frequency. Therefore, 
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changes were detected with the Rodionov method in the SPEI 
(Standardized  Precipitation  Evapotranspiration Index) and  SPI 
(Standardized Precipitation Index) indices at 6 and 12 month scales in 
the state of Durango, Mexico. We analyzed the trends of changes in the 
mean (RSI) and variance (RSSI) of the indices in 49 stations of the SPI 
index for 6 months and in 45 of the SPEI index for 6 months. The RSI 
values identified the beginning of the drought periods, most of the 
maximum RSI values were presented in July 1992, Trends in the 
frequency of RSI values decreased during the 2000-2004 period and 
increased during 2005-2008. The decade of greatest changes in the 
average was 1990-1999. The highest frequency of RSSI values was in 
the 1985-1989 period, the changes have remained constant during the 
last 20 years. The decade of greatest change in variance was 1980- 
1989. 


Keywords: Durango, drought, SPEI, SPI, Rodionov, mean changes, 


changes of variance. 
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Introducción 


La sequía es un fenómeno climático que ha afectado a la humanidad a 
través de la historia (Heim, 2002), con mayor impacto en las personas 
en el planeta (Ortega, 2013). El Instituto Nacional de Ecología y Cambio 
Climático (INECC) indica que la sequía es de los eventos más costosos 


para el país (Ibarrarán €: Rodríguez, 2007). 


Los impactos de la sequía son muy variados, desde la muerte de 
plantas y animales, e inclusive la de seres humanos (Cenapred, 2002), 
efectos sociales en México, con un marcado efecto en el Norte de México 
(Domínguez, 2016). Una sequía no se puede predecir, pero sí es posible 
planear y manejar de modo adecuado los recursos hídricos, a fin de 
mitigar el posible impacto de los daños que pueda causar la sequía 
(Cancelliere € Salas, 2010). 


Así, es necesario un análisis del fenómeno espacial y temporal, 
entre las opciones existen índices de sequía, que tienen por objetivo la 
caracterización de ésta bajo diferentes puntos de vista (Hounam, 
Burgos, Kalik, Palmer, 8 Rodda, 1975). 


Los índices de sequía facilitan el análisis cuantitativo sobre la 


magnitud (con una escala) de la sequía y con ello analizar la distribución 
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espacial y temporal del fenómeno, e incluso su frecuencia. Dichos 
estudios permiten hacer un análisis, evaluación sobre sectores más 
vulnerables, y contar con información útil para la gestión de riesgos en 
los grupos afectados (OMM, 2016). 


En este marco se deben considerar los cambios en el sistema 
climático global que se asocian con condiciones de sequía de mayor 
intensidad, y de mayor duración y frecuencia (Bates, Kundzewicz, Wu, 8 
Palutikof, 2008). Si bien estas variaciones en el clima son ciertamente 
naturales, algunos componentes podrían asociarse con efectos 


antrópicos (Miller, Cayan, € Barnett, 1994). 


En la década de 1970, la interpretación de las variaciones en los 
registros climáticos derivó en la aceptación del término “cambios de 
régimen” (Rodionov, 2005a), que son las reorganizaciones rápidas de 
los ecosistemas de un estado relativamente estable a otro (Rodionov 8 
Overland, 2005). 


Se han desarrollado diversos métodos para detectar cambios de 
régimen, o discontinuidad, en series de tiempo. Típicamente, estos 
métodos emplean técnicas estadísticas estándar, como las pruebas de 
Student o Mann-Kendall, o sus modificaciones (Rodionov, 2004). Un 
problema común de los métodos existentes para la detección de 
cambios de régimen es su bajo rendimiento al final de las series 
temporales. En consecuencia, los cambios en los índices ambientales por 


lo general se detectan mucho después de su aparición real (Rodionov 81 
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Overland, 2005). Como solución a este problema, Rodionov (2004) 
propone el uso de una técnica de procesamiento de datos secuencial, en 
la cual se realiza una comprobación para determinar si una nueva 
observación representa una desviación estadísticamente significativa del 
valor medio del régimen actual y puede considerarse como un cambio 


de régimen. 


En este trabajo se evalúa la sequía a través del Índice de 
Precipitación Estandarizada (SPI), y el Índice de Precipitación y 
Evapotranspiración Estandarizado (SPEI) en el estado de Durango, 
México (con áreas muy afectadas por la sequía). Se aplicó la prueba de 
Rodionov para detectar cambios de régimen de media y varianza de los 
índices de sequía, e identificar si existe tendencia de un incremento o 


disminución de la ocurrencia de cambios en los últimos años. 


Materiales y métodos 


Descripción del área de estudio 
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El estado de Durango se ubica en el noroeste de México y tiene una 
superficie de 123 367 km?, cuenta con 39 municipios (INEGI, 2017). 
Con cuatro provincias fisiográficas (Semarnat, 2009), de mayor a menor 
extensión están la Provincia Sierra Madre Occidental (porción central y 
occidental de la entidad); las Llanuras del Norte (en la zona Noreste); la 


Mesa del Centro (parte centro-oriental), y la Sierra Madre Oriental 
(Figura 1). 


Simbología 

[> División municipal 
[División estatal 
Provincias fisigráficas 


ES LLanura Costera del 
Pacífico 
| Mesa del Centro 


1 Sierra Madre Occidental 

| Sierra Madre Oriental 

- Sierras y Llanuras del 
Norte 


Figura 1. Ubicación del estado de Durango, México, y provincias fisiográficas. 
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Por su localización geográfica, Durango presenta climas semisecos 
templados, secos, muy secos, semifríos, semicálidos y cálidos. La 
distribución de las isoyetas es de suroeste a noroeste en orden 
decreciente, siendo la mayor de 1 200 mm y la de menor de 200 mm; 


esta última ubicada en el límite con el estado de Coahuila (Figura 2). 


Simbología 


[7 División estatal 
Isoyetas de 
—— precipitación 
media anual 
Unidades Climáticas 
2 Cálido subhúmedo 
o Frio 
Muy seco cálido 
Muy seco semicálido 
Seco semicálido 
Seco templado 
- Semifrio subhúmedo 
Semiseco muy cálido 


Sistema de coordenads: 
MEXICO ITRF 2008 LOC Semiseco semicálido 


rojección: 
Lambert Conformal Conie 5 
Datum: ITRF 2008 Semiseco templado 


Templado subhúmedo 


Figura 2. Distribución climática e isoyetas del estado de Durango, México. 


En la vertiente del Pacífico, región que corre de noroeste a 
suroeste, se considera como la zona con mayores precipitaciones (lluvia 
promedio de 1 200 mm), que coincide con las zonas de bosque y selva 


de Durango. En esta parte del estado se encuentra la Sierra Madre 
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Occidental, la cual, al recibir la humedad del Pacífico, y al condesar por 
enfriamiento al ascender, conforma los eventos más importantes de 


lluvia. 


Sólo una pequeña parte de la humedad del Pacífico logra traspasar 
la barrera orográfica hacia la vertiente interior de Durango, en la zona 
de pastizales, matorral y agricultura. Con una precipitación media de 
600 mm anuales, y cada vez menores, conforme se adentra a la parte 
oriental, donde la precipitación media es de 200 mm al año, se 


presentan los climas semiseco, seco y muy seco (Semarnat, 2009). 


Información meteorológica 


La información meteorológica de las variables de temperatura máxima y 
mínima, y precipitación a escala diaria de las estaciones meteorológicas 
que se ubican dentro del área de estudio, se recuperó de las bases de 
datos del Servicio Meteorológico Nacional (SMN). Se descargó la 
información mensual por estación del sitio web 


smn.cna.gob.mx/es/climatologia/informacion-climatologica. 
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Se eligió como periodo de estudio los años de 1979 a 2008, ya 
que incluye un mayor número de estaciones con registros (con al menos 
un 90% y algunas con 75% de registros en zonas con poca información) 
en el periodo de estudio. Con un total de 61 estaciones meteorológicas 
del SMN, de las cuales 52 pertenecen al estado de Durango, tres a 


Coahuila, tres a Sinaloa, dos al estado de Zacatecas y una a Chihuahua 
(Figura 3). 
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Figura 3. Estaciones meteorológicas SMN seleccionadas. 
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A la información meteorológica se le aplicó un control de calidad 
de datos, que consistió en identificar valores no razonables de acuerdo 
con los siguientes criterios: (1) cantidad de precipitación diaria menor 
que cero, y (2) temperatura máxima diaria menor que la temperatura 
mínima diaria. Se eliminaron los valores no razonables y se trataron 


como datos faltantes. 


Los datos faltantes se completaron de dos formas. El primero, con 
el método del US National Weather Service (NWS) —conocido también 
como el cuadrado de la distancia inversa— es un método sencillo de 
utilizar, pero se debe contar con información de estaciones cercanas 
(Campos, 1998); puede ser aplicado a datos diarios, mensuales o 
anuales. Trabajos recientes indican que es uno de los mejores métodos 
para estimar datos faltantes en México (Gallegos, Arteaga, Vázquez, 8 
Juárez, 2016; Gamboa, 2015). El segundo, es el modelo Climate 
Forecast System Reanalysis (CFSR), que desarrolló el National Center 
for Enviromental Prediction (NCEP) y el National Center for Atmospheric 
Research (NCAR) de EUA. Se trata de un modelo que representa la 
interacción global entre los océanos, la tierra y la atmósfera (Saha et 
al., 2010), se cuenta con un periodo de 35 años (1979-2014) de datos 
diarios para las variables de temperatura, precipitación, viento, 
humedad relativa y energía solar; la resolución del CFSR es de unos 38 
km (Saha et al., 2014). 
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El método NWS estima el valor faltante en la estación de acuerdo 
con la Ecuación (1): 


YN Pi*Wi 
Px = —===— 
Y Wi 
(1) 
Donde Pi es la precipitación observada en la estación auxiliar ¡; ¡ = 


1, 2, 3..n, para la fecha faltante; n es el número de estaciones 


auxiliares (mínimo n = 2). 


(2) 


Di es la distancia entre cada estación auxiliar y la estación 


incompleta, en kilómetros. 


Para cada estación se ubicaron las más cercanas como auxiliares y 
se calculó la distancia entre ellas. Cuando fue posible, se ubicaron las 
cuatro estaciones auxiliares más cercanas, una sobre cada cuadrante de 


los ejes coordenados norte-sur y este-oeste (Campos, 1998). 


Los datos del modelo CFSR se obtuvieron del sitio web 


globalweather.tamu.edu. Se realizó la descarga de valores de la malla 
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de puntos para el estado de Durango. Para completar los datos 
faltantes, se eligió el punto de la malla más cercano a cada estación 


meteorológica (Figura 4). 
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Zacatecas 


24D0'N 


Simbología 
LL División municipal 
[7] División estatal 
e Estación meteorológica 


+ _ Punto malla CFSR 
102"b0"w 


106”0'0"W 104*0'0"W 


Figura 4. Puntos de la malla CFSR seleccionados. 


Se completaron los datos faltantes en 59 estaciones con el 
procedimiento para cada variable y el método más adecuado. Para lo 


cual se evaluó la capacidad explicativa de los modelos sobre los datos 
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medidos en las estaciones meteorológicas, a través de cincos índices 
estadísticos: cuadrado medio del error (MSE) (Sánchez-Algarra, 
Barraza-Sánchez, Reverter-Comas, 8 Vegas-Lozano, 2006); raíz 
cuadrada del cuadrado medio del error (RMSE) (Lehmann € Casella, 
1998); error absoluto medio (MAE) (Fuertes 8 Javier, 1998); índice de 
Willmott (d) (Willmott, 1982), y el coeficiente de determinación (R?) 
(Infante € Zarate, 2012). 


Cada índice tiene un criterio para determinar si el modelo es un 
buen estimador; de esta forma se eligió como mejor método para 
estimar datos faltantes aquel en el que al menos tres Índices 


estadísticos lo señalaran como el mejor estimador. 


El modelo CFSR para estimar temperatura máxima en 34 
estaciones resultó ser el mejor estimador de datos perdidos en 
comparación con las 25 que se completaron con el método NWS. El 
método NWS es el mejor estimador para precipitación, pues se evaluó 
como mejor en 57 estaciones, mientras que el modelo CFSR fue el 
mejor sólo en dos. Para temperatura mínima, el método NWS se evaluó 
como mejor en 41 estaciones. De esta manera se obtuvieron 61 series 
climáticas mensuales completas para las variables temperatura máxima 


y mínima, y precipitación. 
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Metodología 


Índices de sequía 


Se utilizaron los índices de sequía Índice de Precipitación Estandarizada 
(SPI) y el Índice de Precipitación y Evapotranspiración Estandarizado 
(SPEI) para el análisis de la sequía. De acuerdo con las 
recomendaciones de la Organización Meteorológica Mundial (OMM, 
2016), ambos ya han sido probados en el área de estudio (Rivera, 
Crespo, Arteaga, 8 Quevedo, 2007; Llanes et a/., 2018) y se pueden 


aplicar a partir de la información disponible. 


El SPI se desarrolló por Mckee, Doesken y Kleist (1993), con el 
objetivo de cuantificar los impactos del déficit de precipitación sobre los 
diferentes recursos hidrológicos: agua subterránea, agua almacenada, 
corrientes de agua y humedad del suelo. El índice se define con base en 
la precipitación estandarizada, que es la diferencia de precipitación de la 
media para un periodo de tiempo específico dividido por la desviación 
estándar, donde la media y la desviación estándar se determinan a 


partir de registros pasados. 
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Se calcula con base en un conjunto de datos mensuales de 
precipitación, el cual se está moviendo en el sentido de que cada mes se 
determina un nuevo valor a partir de los meses anteriores. Cada uno de 
los conjuntos de datos se ajusta a la función Gamma, para definir la 
relación entre probabilidad y precipitación. Dado que la función Gamma 
en análisis climáticos es su integral, se obtienen probabilidades de 
ocurrencia de una precipitación menor o igual que una precipitación 
determinada. Con el propósito de facilitar la obtención del índice SPI, la 
probabilidad acumulada se transforma a una variable normalizada Z 


(con media cero y varianza uno), que representa el valor de SPI. 


El SPEI fue desarrollado por Vicente, Beguería y López (2010). Es 
un índice de sequía multiescalar, que combina datos de precipitación y 
temperatura. Funciona a partir del procedimiento original de cálculo del 
SPI. El SPEI utiliza la diferencia mensual (o semanal) entre la 
precipitación y evapotranspiración (PET). Esto representa un balance 
hídrico climático simple que se fija en diferentes escalas de tiempo para 
obtener el SPEI. 


El cálculo de PET resulta difícil por la cantidad de parámetros que 
intervienen en el proceso. El objetivo de incluir PET en el cálculo del 
índice de sequía es obtener una estimación temporal relativa; por lo 
tanto, el método que se utiliza para calcular el PET no es crítico, es 
decir, que SPEI no se encuentra vinculado con ningún método en 


particular. 
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Para calcular SPEI se necesita una distribución de tres parámetros, 
ya que puede tener valores negativos, que son comunes en la serie de 
balance entre precipitación y evapotranspiración. Vicente et a/. (2010) 
seleccionaron la función de densidad de probabilidad Log-logística, con 
base en el comportamiento de los valores más extremos. Después es 
transformada a una distribución normal estándar con media O y varianza 
1. Dado que SPEI es una variable estandarizada, se puede comparar con 


otros valores de SPEI a lo largo del tiempo y el espacio. 


Mckee et al. (1993) utilizaron un sistema de clasificación para 
definir intensidades de sequía (Tabla 1) y definieron criterios para 
cuando ocurre un evento de sequía, para cualquier escala de tiempo. Ya 
que SPEI es una variable estandarizada, al igual que SPI, se utiliza el 


mismo sistema de clasificación de sequías para ambos índices. 


Tabla 1. Clasificación de SPI y SPEI. 


Valor del índice Categoría del índice 
2 0 mayor Extremadamente húmedo 
1.5a1.99 Muy húmedo 
1.0 a 1.49 Moderadamente húmedo 

-0.99 a 0.99 Cercano a lo normal 
-1.0a-1.49 Moderadamente seco 
-1.5a -1.99 Muy seco 
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Valor del índice Categoría del índice 


-2 0 menor Extremadamente seco 


Se calcularon los índices SPEI y SPI para las series climáticas con 
el programa RStudiof*) (RStudio Inc., 2018), con la paquetería SPEI 
(Begueria, Serrano, 8 Sawasawa, 2013). Se obtuvieron los índices SPEI 
y SPI a escalas de 6 y 12 meses; además de la información climática de 
temperatura máxima y mínima, y precipitación, se utilizó la latitud de 
cada estación. Se seleccionó el método de Hargreaves para estimar la 


evapotranspiración. 


Comprobación de normalidad 


El método de Rodionov utiliza la prueba t-student como parte del 
procedimiento en la detección de cambios. Para emplear el 
procedimiento de la prueba t se debe cumplir el supuesto de que los 


datos proceden de una distribución normal (Montgomery, 2010). 
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Para comparar la normalidad de los datos se eligió la prueba de 
bondad, la prueba de Kolmogorov-Smirnov modificada por Lilliefors 
(1967). Se aplicó la prueba a cada serie de índices de sequía; se realizó 
mediante la paquetería nortest de RStudioO); se utilizó un nivel de 


significancia (alfa, a) de 0.01. 


Método de Rodionov para detectar cambios de régimen 


Rodionov (2004) propuso un método secuencial para detectar cambios 
de régimen (en la media y desviación estándar) como solución al 
problema de las series temporales de no consideran las datos finales de 
los métodos existentes para la detección de cambio de régimen y que en 


consecuencia se detectan mucho después de su aparición real. 


El método consiste en realizar una prueba para determinar la 
validez de la hipótesis nula Ho (en este caso, la existencia de un cambio 
de régimen). Hay tres posibles resultados de la prueba: aceptar Ho, 


rechazar Ho o seguir realizando la prueba. 
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En una serie de tiempo, cuando llega una nueva observación —que 
lo hace de forma regular—, se hace una comprobación para determinar 
si representa una desviación estadísticamente significativa del valor 
medio del régimen actual. Si es el caso, ese año se marca como un 
posible punto de cambio c, y las observaciones posteriores se utilizan 


para confirmar o rechazar esta hipótesis (Rodionov € Overland, 2005). 


La hipótesis se prueba al usar el índice de cambio de régimen 


(RSI), que se calcula para cada c: 


(3) 


Donde m=0,...,1—1 es el número de años desde el inicio de un 
nuevo régimen; 1, la longitud de corte de los regímenes a probar, y o, es 
la desviación estándar promedio para todos los intervalos de un año en 
las series de tiempo. RSI representa una suma acumulativa de 
desviaciones normalizadas x,* del nivel medio hipotético para el nuevo 
régimen (X»vo), para el cual la diferencia del nivel medio para el régimen 
actual (X,..) es estadísticamente significativa de acuerdo con una prueba 
t de Student. 
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dif = Xnvo — Xact = 8 ¡207 /1 


(4) 


Donde t es el valor de la distribución t-student con 2/-2 grados de 
libertad a un nivel de probabilidad p dado. Si en cualquier momento 
desde el comienzo del nuevo régimen el RSI se vuelve negativo, la 
prueba falla y se asigna un valor cero. Si el RSI permanece positivo 
durante /— 1, entonces c se explica como el tiempo de un cambio de 
régimen en el nivel < p. La búsqueda del próximo cambio de régimen 
comienza con c+ 1 para garantizar que su sincronización se detecte de 


forma correcta, incluso si la duración real del nuevo régimen es < 1 año. 


El valor promedio del régimen actual x,.. se calcula para el periodo 
(c— lc). Se calcula para el periodo desde el cambio de régimen anterior 
al punto inmediatamente anterior al punto actual en el tiempo. Como 
resultado, se produce una función gradual de los regímenes en casi 


todos los casos y no sólo para cambios abruptos en el régimen. 


Los cambios en la varianza de alguna variable meteorológica 
pueden tener un impacto similar o incluso mayor en los ecosistemas que 
las variables en la media. Como los cambios climáticos debidos a causas 


naturales o humanas pueden ser en la frecuencia de los eventos 
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extremos, representan el peligro de favorecer el aumento del 


"calentamiento global" (Rodionov, 2005b). 


Rodionov desarrolló una extensión para Excel (Rodionov, 
2005b) que detecta cambios en la varianza, procedimiento similar al de 
la media, excepto que es con base en la prueba F en lugar de la prueba 
E. 


Los datos de entrada para el programa son los valores de los 
índices de sequía y los parámetros de la prueba: longitud de corte, nivel 


de significancia y función de peso de Huber. 


Los parámetros de longitud de corte y nivel de significancia 
controlan la magnitud y escala de los regímenes que fueron detectados. 
Diferentes combinaciones en estos parámetros tienen un efecto directo 
sobre el valor de RSI, es decir, se obtendrán diferentes interpretaciones 
de los resultados (Rodionov, 2006a). De esta manera, se analizó el caso 
de una longitud de corte igual a cinco para un nivel de significancia de 
0.05. 


El parámetro de peso de Huber controla la participación de los 
valores atípicos. Debido a estos valores, el promedio no es 
representativo del valor medio de los regímenes. El peso para el valor 
de los datos se debe elegir de modo que sea pequeño si ese valor se 
considera como un valor atípico (Rodionov, 2006a). En función de la 


desviación estándar se usó un peso de Huber de 2. 
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Análisis de tendencias de cambios en media y varianza 


de los índices de sequía 


Para cada estación se calcularon las frecuencias de ocurrencia de 
cambios en la media (RSI) y varianza (RSSI) de los índices de sequía; 
para los valores RSI se calculó sólo la frecuencia de los valores RSI 


negativos, que representan los regímenes de sequía. 


Se hizo un análisis espacial de las frecuencias obtenidas de los 
valores RSI y RSSI a través de la interpolación, que fue con base en el 
principio de autocorrelación y determina si hay un patrón espacial 
(Childs, 2004). La interpolación se hizo con la extensión Spatial Analyst 
de ArcGISO (ESRI, 2016), con el método kriging (Villatoro, Henríquez, 
8: Sancho, 2008). 


Para el análisis de tendencias, las frecuencias de ocurrencia de 
cambios en la media y varianza se agruparon en periodos de 5 y 10 
años. Se identificó si durante estos periodos se presentó una tendencia 


de disminución o incremento del número de cambios. 
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Resultados 


Se calcularon los índices de sequía SPEI y SPI de seis meses para el 
periodo de junio de 1979 a diciembre de 2008; para los índices SPEI y 
SPI de seis meses, los valores obtenidos son de diciembre 1979 a 
diciembre de 2008. En la Tabla 2 se presentan las estadísticas de los 


índices de sequía. 


Tabla 2. Estadísticas generales de los índices de sequía para el estado de Durango, 


México. 
P Desviación , ] 
Indice de sequía | Media E Mínimo Máximo | Rango 
estándar 
SPI de 6 meses 0.009 0.984 -3.025 2.46 5.485 
SPI de 12 meses 0.003 0.986 -2.752 2.376 5.128 
SPEI de 6 meses 0.008 0.979 -2,308 2.268 4.576 
SPEI de 12 meses 0.006 0.978 -2.157 2.161 4.318 
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La media de los índices (SPI y SPEI) es muy cercana a cero, lo 
cual indica que los valores extremos (periodos de sequías y de 
humedad) ocurren con la misma frecuencia. Es decir, que la presencia 
de las sequías durante el periodo de estudio es una característica normal 
del clima de la región, resultado similar a lo que presentó Contreras 
(2005). 


En términos de la magnitud de los promedios de valores extremos, 
es mayor para los valores mínimos en tres de los índices analizados (SPI 
de 6 y 12 meses, y SPEI de seis meses). Esto indica que la magnitud de 
la sequía es mayor con respecto a la magnitud de los periodos húmedos. 
Ello se traduce en que las sequías representan mayor riesgo en el área 


cuando se presentan. 


De las pruebas de normalidad, se cumplieron en 46 estaciones 
para el índice SPI de seis meses; 45, para SPEI de seis meses; 32, del 
índice SPI de 12 meses, y 31 para SPEI de 12 meses. Se realizó la 
prueba de Rodionov para la detección de cambios de régimen para los 
índices SPI y SPEI de seis meses, pues son los que cuentan con un 
mayor número de estaciones que cumplen la condición de normalidad, 


esto, con la finalidad de incluir una mayor superficie en el análisis. 
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Análisis del cambio de régimen en la media (RSI) y la 
varianza (RSSI) 


Para la media de los índices SPI y SPEI de seis meses se obtuvo el RSI; 
la media de los regímenes con pesos iguales; el valor medio de los 
regímenes calculados con la función de ponderación de Huber; la 
longitud de los regímenes; el nivel de confianza final para los cambios, y 


el peso (en función de la desviación estándar) de los valores atípicos. 


Para la varianza, se obtuvo el índice de cambio de régimen de la 
varianza promedio (RSSI) y, al igual que para la media, se obtuvo la 


longitud del régimen y el nivel de confianza para los cambios. 


Con los valores de RSI y RSSI se calcularon las frecuencias de 
ocurrencia de los valores de cambio por estación para el periodo 1979- 
2008. 


Para RSI se obtuvieron las frecuencias de los valores negativos, 
que son los que indican el cambio hacia un régimen de sequía. Para 
cada estación se obtuvieron los valores mínimos y máximos de RSI y 


RSSI, respectivamente, así como la fecha en que se presentaron. 


Los valores mínimos promedio de RSI son -3.23 para el índice 


SPEI de seis meses y -2.96 para SPI de seis meses. El promedio de 
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valores máximos de RSSI es de 1.30 para SPEI de seis meses y 2.00 


para SPI de seis meses. 


Para los cambios en la media se tiene mayor frecuencia de 
ocurrencia de valores máximos en una fecha específica; mientras que 


para los cambios en la varianza, la frecuencia es más constante. 


A través del análisis espacial se identificaron áreas del estado que 
muestran mayor frecuencia de cambios de régimen en la media y 
varianza. Para el RSI, en ambos índices, los valores de menor ocurrencia 
se presentan al sur del estado, que comprende un área mayor en el 
índice SPI. Una estación que se ubica en el noroeste presentó la mayor 
frecuencia para el índice SPEI; mientras que para el índice SPI, la 


frecuencia de cambios RSI fue la menor (Figura 5). 
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Figura 5. Distribución espacial de la frecuencia de valores RSI durante el periodo 


1979-2008 en el estado de Durango México. 


En la Figura 6 se observa que la frecuencia de los valores de 
cambio en la varianza RSSI va de 5 a 13 durante el periodo 1979-2008. 
Para el índice SPEI de seis meses, la menor frecuencia se presenta al 
oeste del estado, a modo de una franja que va de norte a sur, 
excluyendo una pequeña porción en los límites con Chihuahua y Sinaloa, 
que presenta una frecuencia de 7 a 9 ocurrencias, al igual que la mayor 
parte del estado, ya que las frecuencias mayores sólo representan una 
correspondiente con valores de sólo dos 


pequeña área al este, 


estaciones climatológicas. 
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Figura 6. Distribución espacial de la frecuencia de valores RSSI durante el periodo 


1979-2008 en el estado de Durango, México. 


La frecuencia de RSSI para el índice SPI de seis meses presentó 
dos categorías de ocurrencia de 7 a 9 que cubre la mayor parte del 
estado; de 9 a 11 distribuida al oeste (resultado de dos estaciones), y 


una pequeña porción al este (de una estación). 


La detección de cambios en la media y varianza depende del índice 
de sequía utilizado. Al emplear el índice SPEI, las frecuencias muestran 
mayor variabilidad, al agruparse en mayor número de categorías que las 


frecuencias obtenidas para el índice SPI. 


Tecnología y ciencias del agua, ISSN 2007-2422, 11(1), 85-131. DOI: 10.24850/j-tyca-2020-01-03 


114 


2020, Instituto Mexicano de Tecnología del Agua 
Tecnología y 


ienciaszAgua Open Access bajo la licencia CC BY-NC-SA 4.0 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


De acuerdo con el análisis espacial, la región norte del estado 
presenta mayor frecuencia de cambios; la estación que muestra mayor 


número de cambios es la 10078 Sardinas. 


Por lo cual, a continuación se presenta un análisis individual de los 
resultados obtenidos en esta estación para los valores de RSI y RSSI 
(cambios en media y varianza de los índices de sequía, 


respectivamente). 


Se representa en la gráficas (Figura 7) los valores de los índices 
de sequía, así como los valores de RSI y RSSI. Para los valores de RSI 
negativos, se ubicaron en el inicio de los periodos de sequía para los 
índices SPEI y SPI de seis meses, es decir que, el método detectó todos 


los cambios hacia un régimen de sequía. 
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Estación 10078 Sardinas. Índice SPEl 6 meses, 1979-2008 
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Figura 7. Valores de cambio RSI para los índices SPI y SPEI de seis meses en la 


estación10078 Sardinas. 


Hubo 18 números de cambios en la media para el índice SPEI de 
seis meses, 18 en total; y 16 cambios para SPI en el mismo periodo. La 
mayor parte de los valores RSI se calculó en la misma fecha para ambos 


índices. La detección de cambios en la media depende entonces del 
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índice de sequía con el cual se trabaje (SPEI o SPI), siendo más sensible 


con el índice SPEI. 


De 1994 a 1996 se presentó el mayor periodo de sequía en la 
estación 10078 Sardinas. En este periodo se observan dos cambios 
negativos en la media. El valor de RSI es útil para detectar si una sequía 
continuará, pues el valor negativo de RSI indica que los valores de 


clasificación de sequía de los índices seguirán en la misma tendencia. 


Los valores mayores de RSI se tienen cuando el cambio de 
régimen es más abrupto; en cambio, cuando la media es más constante, 
el cambio se da de forma gradual; los valores de RSI son de menor 
magnitud. Si los índices de sequía se encuentran en valores grandes de 
humedad (positivo) o de sequía (negativo), el cambio a un valor grande 


de la condición contraria resulta en valores de RSI de mayor magnitud. 


Los valores máximos de RSI ocurrieron en diversas fechas en las 
estaciones analizadas, sin embargo se detectó que una gran parte 


ocurrió durante los meses de junio y julio de 1992. 


Los cambios de la varianza no indican un comportamiento general 
de la sequía, sino que es particular de acuerdo con el índice utilizado. 
Además, la varianza es sensible a los valores extremos, por lo que la 
interpretación de los resultados se inclina hacia un cambio rápido en los 


valores de los índices de sequía. 
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Los cambios detectados en la varianza para la estación 10078 
Sardinas (Figura 8) indican una varianza media constante en 
prácticamente todo el periodo. Para SPEI de seis meses, los cambios en 
la varianza se ubican en los años 1985 a 1987; para SPI, la mayor parte 


de las variaciones se detectaron al final de la serie, de 2005 a 2007. 


Estación 10078 Sardinas. Índice SPEl 6 meses, 1979-2008 
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Figura 8. Valores de cambio RSSI para los índices SPI y SPEI de seis meses en la 


estación 10078 Sardinas. 


En la estación 10078 Sardinas, los valores máximos de RSSI no 
coinciden con los valores máximos de RSI; en cambio, los máximos de 
RSSI se observan en una variación de la media de los índices de sequía, 
lo cual es de esperarse, ya que el índice de varianza indica la dispersión 


de los valores con respecto a su media. 


Con lo anterior se puede inferir que los valores de RSSI serán 
útiles para detectar cambios hacia un periodo prolongado de sequía. Sin 
embargo, cuando los periodos de humedad se alternen con periodos de 
sequía con la misma magnitud y duración, no se detectarán cambios en 
la varianza. Por lo cual, los valores de RSSI se deben utilizar como un 
complemento de los valores RSI. Con ambas pruebas es posible que se 
detecten con mayor certeza los periodos de sequía prolongados y de 


mayor severidad y duración. 


Se analizó la frecuencia de los valores de cambio RSI y RSSI de la 
media y varianza en los índices de sequía mediante la suma de valores 
cada 5 y 10 años. De esta manera, se identificó la tendencia de la 
frecuencia de cambios; es decir, si ha aumentado o disminuido la 
ocurrencia de cambios en la media y varianza de los índices de sequía 


en el tiempo. 
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La tendencia de la frecuencia de cambios en la media cada cinco 
años (quinquenal) para los índices de sequía SPEI y SPI de seis meses 
se indica en la Figura 9. La frecuencia de ambos índices es parecida; 
sólo durante el periodo 1995-1999 la frecuencia de cambios para el 
índice SPEI de seis meses es menor que para SPI de seis meses. 
Durante los primeros cuatro quinquenios se observa una tendencia de 
aumento en los cambios en la media, y después durante el periodo 
2000-2004 disminuye radicalmente; en el último periodo de 2004-2008, 
los cambios aumentan de nuevo hasta valores similares que en los 


primeros quinquenios. 
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Figura 9. Frecuencia quinquenal de cambios en la media (RSI) durante el periodo 
1980-2008. 
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La frecuencias de cambios en la media cada diez años (decadal) se 
tiene en la Figura 10. Se observa que durante la última década (2000- 
2008) ocurrió la menor cantidad de cambios, y durante el periodo 1990- 


1999 se presentó la mayor frecuencia de valores RSI. 
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Figura 10. Frecuencia decadal de cambios en la media (RSI) durante el periodo 1980- 
2008. 


Los resultados obtenidos coinciden con los encontrados por Llanes 


et al. (2018), donde se concluye que los eventos extremos de sequía se 
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presentaron durante los años de 1990 y 1997, siendo los años 


intermedios un periodo de sequías recurrentes. 


Para la frecuencia de cambios en la varianza, los índices de sequía 
SPI y SPEI de seis meses se muestran en la Figura 11, donde se 
observa una mayor frecuencia (> 40) de cambios para el índice SPI de 
seis meses en todos los periodos; la frecuencia de RSSI es mayor para 
el SPI con respecto a SPEI. 
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Figura 11. Frecuencia quinquenal de cambios en la varianza (RSSI) durante el periodo 
1980-2008. 


La frecuencia de cambios en la varianza para periodos de diez 


años se muestra en la Figura 12. Se observa que la primera década de 


Tecnología y ciencias del agua, ISSN 2007-2422, 11(1), 85-131. DOI: 10.24850/j-tyca-2020-01-03 


122 


2020, Instituto Mexicano de Tecnología del Agua 
Tecnología y 


CienciaszAgua Open Access bajo la licencia CC BY-NC-SA 4.0 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


análisis (1980-1989) es la que presentó mayor frecuencia de cambios; 
sin embargo, para SPI, el valor se repite en la siguiente década, por lo 
que durante 20 años el promedio de cambios fue relativamente 
constante. 
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Figura 12. Frecuencia decadal de cambios en la varianza (RSSI) durante el periodo 
1980-2008. 


Los resultados obtenidos indican que los cambios en media y 
varianza no han sido mayores durante los últimos años del periodo 


analizado. 
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Conclusiones 


Se analizó la sequía por medio de los índices de sequía SPI y SPEI a 
escalas temporales de 6 y 12 meses en el estado de Durango, México. 
Se identificaron los periodos de sequía con los valores de cambios de la 
media (RSI) de signo negativo. Se detectaron 14 cambios en promedio 
para el índice SPEI de seis meses y 13 cambios para el índice SPI de 
seis meses. Los valores de cambio en la varianza (RSSI) fueron ocho 


para el índice SPEI de seis meses y 10 para el índice SPI de seis meses. 


Las tendencias de los cambios en la media y varianza de los 
índices de sequía con base en la frecuencia de los cambios en la media 
disminuyó durante el periodo 2000-2004 y aumentó durante 2005- 
2008. La década de mayores cambios en la media fue 1990-1999. Para 
la varianza se observó una frecuencia máxima de cambios durante el 
periodo 1985-1989; los cambios se mantuvieron constantes durante los 
últimos 20 años. La década de mayores cambios en la varianza fue de 
1980-1989. 
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